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Photothernial raicroscopy [1-4] has recently comple- 
mented single molecule fluorescence microscopy by the 
detection of individual nano-objects in absorption. Pho- 
tothermal techniques gain their superior sensitivity by 
exploiting a heat induced refractive index change around 
the absorbing nano-object. Numerous new applications 
to nanoparticles, nanorods and even single molecules [5- 
10] have been reported all refering to the fact that pho- 
tothermal microscopy is an extinction measurement on 
a heat induced refractive index profile. Here, we show 
that the actual physical mechanism generating a pho- 
tothermal signal from a single molecule/particle is fun- 
damentally different from the assumed extinction mea- 
surement [11, 12]. Combining photothermal microscopy, 
light scattering microscopy as well as accurate MiE scat- 
tering calculations to single gold nanoparticles, we reveal 
that the detection mechanism is quantitatively explained 
by a nanolensing effect of the long range refractive index 
profile. Our results lay the foundation for future devel- 
opments and quantitative applications of single molecule 
absorption microscopy. 

Long range interactions play a very special role in 
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Fig. 1. Single gold-nanoparticle light scattering, (left) 
Schematic representation of the experimental setup used for 
the images (sample scanning). Shown is the particle posi- 
tion Zp^ whereas Xp and yp are the lateral directions (not 
depicted), (center) Scattering of the heating laser on the gold 
nanoparticle (AuNP). (right) Scattering of the probe laser 
on the AuNP [R — 30 nm). The top-row images represent 
the scans for a detection aperture NAcz = 0.8, the bottom 
row images NA^ = 0.3. The images are grouped experimen- 
tal/theoretical scans. 



physics. Fundamental forces such as the Coulomb force 
or gravitation are attributed to interaction potentials 
with inverse distance dependencies manifesting their ac- 
tion even at macroscopic length scales. In this paper we 
explore the action of another long-ranged field interacting 
with photons, namely the refractive index profile gener- 
ated in the solvent surrounding a suspended nanoparticle 
that releases heat. A typical practical situation is pro- 
vided by a gold nanoparticle in the focus of a laser beam. 
Most of the absorbed optical energy is released as heat 
to the solvent. This creates a temperature field which 
decays with the inverse distance r from the absorber ac- 
cording to 



AT(r) = ATc 



R 



(1) 



where the temperature rise at the particle surface ATq = 
Pahs/^'TThzR is resulting from the absorbed optical power 
Pabs^ the surroundings' heat conductivity k and the par- 
ticle radius R. A corresponding refractive-index profile 
is thereby established, modifying the unperturbed refrac- 
tive index rim of the solvent by (eqn. 2): 
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This infinite hot lens can be exploited to detect even 
a minute absorber with a probe laser in an optical mi- 
croscopy setup with extremely high sensitivity [4]. How- 
ever, the infinite size of the lens and the complex spatial 
structure of the tightly focused, aberrated laser beams 
cause some conceptual and computational challenges for 
a quantitative understanding of the signal generation [12]. 
These are addressed below, where we develop a consis- 
tent mathematical formalism to quantitatively explain 
and analyze such "absorption microscopy" experiments 
in great detail. 

The complexity of the signal generation can already 
be appreciated when both involved lasers (heating and 
probe) scatter from a single gold nanoparticle at low inci- 
dent power, where the heating is still negligibly small (see 
methods section for experimental details). The trans- 
mitted intensity collected by a lens reveals strong inter- 
ference patterns as it senses the field structure of the 
aberrated incident laser beams (Fig. 1, right). Recording 
these scattering intensities at different detection numer- 
ical apertures emphasizes the importance of the phase 
relation of scattered and transmitted electric fields [13]. 
The interference patterns largely vanish when consider- 
ing the photothermal signal, which is the difference in the 
probe laser scattering signal of a heated nanoparticle in- 
cluding the long range refractive index profile and a non- 
heated particle without the refractive index change. This 
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Fig. 2. Photothermal signal as a function of heat- 
ing and probe laser focus displacement. Difference of 
hot/cold NP scattering images: a) Static difference image, b) 
Dynamic PT signal recorded with a lock-in amplifier, c) Com- 
parison of axial PT dynamic (solid black) / static (red dots) / 
theoretical (solid red) scans, d) Schematic of the defocussing 
parameter Azf and the axial particle position Zp. e) Experi- 
mental PT signal traces for different defocussing parameters 

Azf. 



scattering difference has been evaluated under steady 
state conditions recording two scattering images as well 
as dynamically using the photothermal heterodyne tech- 
nique [12] for an axial focus displacement of both lasers of 
350 nm (see Fig. 2 a and b). Despite the much lower sig- 
nal to noise ratio of the steady state difference, both sig- 
nals agree perfectly and reveal a two-lobe structure (Fig. 
2 c). The simplicity of this two- lobe structure suggests 
a much simpler mechanism than the complex aberrated 
spatial heating and probe laser scattering intensity dis- 
tributions put forward. It also unveils that the common 
assumption of a product point spread function of heat- 
ing and detection laser is not appropriate. Instead, the 
two-lobe structure depends sensitively on the displace- 
ment of the two laser foci. Whenever the probe laser 
focus is in front of the refractive index gradient, the de- 
tected intensity is decreased, while an increased signal is 
measured when the probe laser focus is behind the re- 
fractive index gradient, as indicated in Figure 2 e. This 
is exactly the action of a diverging lens. While it is in 
general tempting to consider the refractive index profile 
as a finite sized scatterer which has to be treated by an 
appropriate scattering theory, we infer that the missing 
characteristic length scale of refractive index profile is 
the source of this intuitive lens-like action. As most ma- 
terials lower their refractive index with increasing tem- 
perature (dn/dT < 0), the temperature field induces a 
diverging lens with radial symmetry. In fact a more de- 
tailed analysis reveals that the problem is equivalent to 
Rutherford scattering of a-particles on a Coulomb 
potential in which case a classical mechanics treatment 



of the Qf-particle delivers the same result as a quantum 
mechanical treatment of the particle [26]. While we will 
detail this fundamental equivalence in a separate paper 
[14], we will stress a complex scattering description to 
demonstrate that an exact electromagnetic treatment is 
providing the same picture of a gradient index nanolens. 

This theoretical scattering description has to go be- 
yond common MiE theory as photothermal microscopy 
employs highly focused laser beams instead of plane 
waves. Thus a more rigorous MiE description (Gener- 
alized LoRENZ-MiE Theory, GLMT [24]) has been ex- 
tended to account for the axial structure of the signal 
and to accurately model aberration effects in the focused 
laser beams [15, 16] interacting with the metal particle 
and the refractive index gradient. The latter one is intro- 
duced in a multiple shell approach discretizing the refrac- 
tive index profile as described by Peha et al. [17] (see Fig. 
3 a). The detection aperture NA^^ has been introduced 
into the GLMT formalism calculating the time-averaged 
Poynting vector (S)^ integrated over an area Ad repre- 
senting an angular detection domain which is determined 
by the detection aperture (eqn. 3). 



L 



{S(r))t • dA = - / 7^{E(r) xH*(r)}- dA 

(3) 

Incorporating all experimental and material parameters 
(see method and supplement) to evaluate the pure gold 
nanoparticle scattering, an unprecedented quantitative 
agreement with the experimental results is unveiled (see 
Fig. 1). To our knowledge, this is the first quantitative 
exact calculation of single gold nanoparticle scattering 
intensity distributions. The calculated scattering images 
are largely determined by aberrations of the incident laser 
fields and additional interference structures. The maxi- 
mum detected scattering intensity and true focus position 
do not coincide and strongly depend on wavelength. The 
spatially extended intensity distribution causes the peak 
intensity to be much lower than for a Gaussian beam hav- 
ing a beam-waist of the corresponding diffraction limit 
(^ 0.61A/NAiii). Thus, as compared to such a Gaus- 
sian beam only one third of the particle temperature is 
reached. The amplitude of the aberration induced inter- 
ference pattern is determined by the numerical aperture 
of the detection lens as found in the experiment. The 
larger this detection aperture, the weaker the interfer- 
ence patterns get as the individual phase differences at 
different detection angles average out (Fig. 1 right). This 
finding confirms that the plane-wave optical theorem is 
not applicable when collecting signals at finite detection 
angle and focused illumination [18]. The only way to con- 
clude on the intensity distribution details in the scatter- 
ing experiment is thus to consider the finite detection an- 
gle. While this complex structure determines the shape 
of the scattering intensity distribution, a contribution to 
the photothermal signal is also expected [19, 20] but less 
obvious. 

We evaluate the photothermal signal as the difference 
between the probe scattering intensity distribution of 
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Fig. 3. Generalized Mie scattering calculation results 
of the photot hernial signal, a) Illustration of the temper- 
ature profile T(r), the refractive index profile n (r) and the 
discretization of the latter, b) Theoretical photothermal Zp 
scans for varying foci displacements Azf. c) Peak amplitudes 
vs. Azf. d) Peak positions vs. Azf. e) Signal decomposition 
for do R — 10 nm AuNP at the positive lobe Zp position: Pext 
(black), Psca (blue), Pext + Psca (red). The inset shows Pext 
and Pext + Psca for n (r) = rim + /^nR^ /r'^ (dashed / orange) 
and n (r) — Um + Anexp(— (r — R) /2R) (dashed / red) 



particle with the refractive index profile and the distri- 
bution without the profile, completely equivalent to the 
above presented experiments. Accordingly, the theoreti- 
cal model directly confirms the experimental observations 
quantitatively. Almost all of the interference structure 
disappears due to the difference of the two aberrated sig- 
nals. A two lobe structure is remaining, which exactly 
corresponds to the experimental observations (Fig. 2 c). 
Even the strong dependence of the total signal on the 
displacement of the two involved laser foci is reproduced 
validating our theoretical approach (Fig. 3 b). Thus the 
exact treatment of the signal by a generalized MiE scat- 
tering calculation predicts a lens-like action of the re- 
fractive index gradient as well. It demonstrates that a 
quantitative analysis of photothermal microscopy data 
in terms of temperatures, absorption cross sections and 
even in terms of the sign of the t her morefr active coef- 
ficient dn/dT is possible. This puts photothermal mi- 
croscopy to a new quantitative level. 

The fundamental difference between the commonly as- 
sumed extinction process and the mechanism reported 
here becomes obvious when separating the individual 
contributions to the signal as done for instance by Gaiduk 



[11]. The interaction of a probing field Ep^ with the heat 
induced scatterer is commonly solved through the intro- 
duction of an outgoing spherical wave Egca- The detected 
power Pd then mathematically decomposes, upon insert- 
ing E = Epr + Egca (and H) into eqn. 3, into three sep- 
arate integrals [23] for probe background Ppr, scattering 
Psca and extinction Pext powers containing only the prob- 
ing field Epr, the scattered field Egca and both in mixed 
terms, respectively. So far, it has been commonly as- 
sumed that solely the latter extinction part contributes 
to the photothermal signal [12]. This assumption does 
not hold, independent of the actual gold nanoparticle size 
thus being valid for single molecule detection as well. A 
separation of scattering and interference contribution in 
our theory shows that both parts change in a nonlinear 
fashion when increasing the temperature of the nanopar- 
ticle (see Fig. 3 e). While up to a temperature rise of 
about 50 K the interference contribution determines the 
signal, the importance of the scattering contribution rises 
quickly. It is the consequence of the inverse distance de- 
pendence of the refractive index change, which removes 
the characteristic length scale from the profile and thus 
requires a detailed theoretical analysis. This actually 
becomes obvious when considering a gold particle sur- 
rounded by an artificial refractive index profile which de- 
cays exponentially or with the inverse distance squared. 
In both cases, the finite length scale of such a refractive 
index profile ensures a leading contribution coming from 
the interference term (Fig. 3 e inset) increasing linearly 
with the temperature rise of the particle. Therefore, the 
missing length scale of the refractive index change is re- 
sponsible for both scattering and interference contribu- 
tions to the signal, which have to be treated explicitly to 
understand the signal generation mechanism. The overall 
quality of our theoretical and experimental results is best 
seen in Figure 3 c and d, plotting the maximum and mini- 
mum photothermal signal values as well as their positions 
as a function of the axial focus displacement Azj. They 
show an unprecedented agreement of experiment (mark- 
ers) and theory (lines) without the inclusion of fitting 
parameters. From the above consideration of different 
refractive index profiles it is evident, that a calibration 
of the photothermal signal for the measurement of ab- 
sorption cross-sections on arbitrarily shaped objects of 
sizes comparable to the lateral focus-extent (e.g. single- 
walled carbon-nanotubes [5, 6]) will fail as the general 
contribution of scattering and interference will change. 
In the case of point-like absorbers, however, there is a 
quantitative simplification of the generalized MiE theory 
applied here. In fact the lens-like action can be treated 
in a geometrical optics framework due to its fundamental 
mathematical equivalence to Rutherford scattering as 
we will show elsewhere [14]. 

Finally, the contribution of aberrations to the pho- 
tothermal signal can be rationalized. The photothermal 
signal generated by a non-aberrated probe beam is van- 
ishing if the particle is in the focus of the probe laser. 
This behavior is expected for a simple lens. The aber- 



250 


Az/ = -400nm /^^i 




200 






1 50 






100 






50 













-2-10 1 


2 




Zp [|jm] 





6x10 

4 ^ 

H 

2 W 

CQ 





-3 





Fig. 4- Comparison photothermal and fluorescence 
signal. Left column: Calculated scans of the absorbed power 
^abs (oc point spread function, |E|^) for the heating laser 
(top) and the detection laser (bottom). Right column: Rep- 
resentative scenarios for the relative alignments of the two 
aberrated lasers: (top) maximal positive photothermal sig- 
nal, (center) symmetric signal and (bottom) maximal nega- 
tive signal. Plotted are the heating beam (green, fluorescence) 
and probe beam (red) intensities as well as the resulting pho- 
tothermal signals (blue) vs. the axial particle coordinate Zp. 



rated probe beam, however, leads to a finite signal at the 
particle position. This can be understood in terms of a 
lens-like action as well. While a non-aberrated intensity 
distribution shows a zero crossing of the signal at the 
particle position (not shown), the aberrated probe beam 
leads to a finite signal at the particle position. While a 
non-aberrated beam is symmetric to the lens if directly 
focused to the center of the lens, the aberrated beam is 
not (Fig. 4, left). The additional interference maxima 
act like additional foci displaced with respect to the pho- 
tothermal lens position and thus cause a photothermal 
signal. As a consequence, the displacement of particle 
fluorescence excited by the heating laser and photother- 
mal signal will depend on the aberration of the probe 
beam and not simply on the Gouy phase [21]. We have 
calculated three limiting cases, where either the positive 
lobe or the negative lobe is maximum or both lobes show 
the same photothermal signal magnitude (Fig. 4, right). 
In the first case, photothermal signal maximum and flu- 
orescence signal maximum would be displaced by about 
100 nm (top). In the second case, the fluorescence sig- 
nal is almost exactly at the position of the photothermal 
signal minimum (center) while in the last case the zero- 
crossing of the photothermal signal coincides with the flu- 
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Fig. 5. Principle scheme of the experimental setup. PD: 
photo-diode, P: pinhole, D: dichroic mirror, F: filter, AOM: 
accusto-optic modulator, a: variable ND filter, ADC: Adwin 
analog digital converter. 



orescence maximum (bottom). Therefore the commonly 
found displacements of fluorescence and photothermal 
signal [5, 22] in axial direction can be well explained by 
the presented theoretical approach. 

In summary, our results demonstrate that the long 
range refractive index change generated by a single 
heated nano-object in photothermal microscopy act as 
a nano lens. The understanding of this signal genera- 
tion mechanism establishes photothermal microscopy as 
a quantitative technique to determine absolute absorp- 
tion cross sections and delivers a framework for new ap- 
plications of this technique. As a direct consequence, new 
experimental techniques such as twin-focus photother- 
mal correlation spectroscopy similar to the well estab- 
lished dual focus fluorescence correlation techniques [27] 
or super-resolution absorption microscopy methods can 
be developed. We expect that the understanding of the 
photothermal signal generation will pave the way for fur- 
ther improvements of this technique beyond the current 
level of sensitivity. 



A. Sample preparation. 

Samples were prepared by spin-coating a polymer layer 
(Sylgard 184, about 15 fim thickness) on top of a glass- 
cover slide. AuNPs (BBI International) with a diameter 
of 60 nm were deposited on the polymer film. The parti- 
cles were covered with a second Sylgard layer (also about 
15 /im thick) to embed the particles in a homogeneous 
matrix. The absence of a close-by interface ensures a ra- 
dially symmetric temperature field around the particle. 



B. Photothermal microscopy measurements. 

The experimental setup for single particle light scatter- 
ing and photothermal measurements is based on a home- 
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built confocal sample-scanning microscopy setup using 
two laser sources. A DPSS laser (Coherent, Verdi) with 
Xh = 532 nm is used to heat the gold particles and a sec- 
ond laser source at Xd = 635 nm (Coherent ULN laser 
diode) probes the local refractive index changes. Both 
beams are focused into the sample by the same objec- 
tive lens (Olympus 100x/1.4NA) and are collected above 
the sample by a second objective (Zeiss 10x/0.3NA or 
Olympus 50x/0.8NA), which is adjusted to image the 
probe focus to infinity. The sample is moved by a piezo- 
scanner (PI) . The resulting parallel beam is focused onto 
two photodiodes (Thorlabs, PDA36A-EC) after passing 
appropriate filters (no pinhole). To allow for a low noise 
detection of the photothermal signal, the heating beam 
is modulated with a frequency ft = 300 kHz and the 
probe signal change is detected at this reference fre- 
quency with a lock-in amplifier (Signal Recovery 7280 
DSP) with a time constant of ru = 1 ms. The resulting 
lock-in signal is recorded by a A/D converter (Adwin- 
Gold, Jager Messtechnik) 300 times for each recorded 
pixel (Ims/pixel). The relative photothermal signal ^ 
corresponds to the relative modulation amplitude at the 
photodiode 
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C. Single particle light scattering measurements. 

Single particle light scattering has been carried out 
in the same setup as the photothermal measurements. 
To do so, the intensity of heating and probe laser have 
been diminished that no notable effects of varying inci- 
dent laser power can be found in the experiments. The 
laser intensities of both wavelength were recorded in- 
dependently with two photodiodes (Thorlabs, PDA36A- 
EC) without intensity modulation and lock-in detection. 
All scattering signals have been normalized to the back- 
ground intensity. 



D. Mie scattering calculations. 

Mie scattering calculations were carried out with a 
modified C-code. The multishell scatter coefficients a^+^ 
and b^^^ were provided by ref. [17], and are referred to 
in the supplemented material as dji and bn- The modifi- 
cations include the numerical evaluation of beam shape 
coefficients (supplement) as well as the integration of 



the PoYNTlNG-vector across a finite azimuthal collection- 
angle range (supplement eqn. 5-8). The C-code was then 
interfaced to WaveMetrics Igor Pro 6.06 to create the 
scans and images. 

E. Parameters used for the calculations in Figure 1 

The following parameters were used for the calculation 
of the imaging (subscript d and h denote detection and 
heating parameters, see supplement for notation): Parti- 
cle parameters: R = 30 nm, riAui^h) = 0.516 -h 2. 23i and 
^Au(^d) = 0.175 + 3.46i [25]. Beam and objective param- 
eters: Detection laser: PpM,d = 250/iW, Xd = 635 nm, 
heating laser: PpM,h = 25/iW, Xh = 532 nm, effective il- 
lumination objective focal length / = 1.8mm, overfilling 
factors f/uJd = 7d = 1.125 and f /uj^ = 7/^ = 1. Sam- 
ple parameters: refractive indices: no = 1, ni^h = 1.525, 
m^d = 1.505, nl = 1.514, = = 1.46,' dn/dT = 
-3.6 X IQ-^K-^, hi = 0.15Wm"^K-i, kod"" = 720, 
d = 15 /im, NAiii = 1.4, NA^ = 0.3 or NA^ = 0.75. Dis- 
cretization parameters: number of angles for numerical 
angular signal integration: Nq = 400, number of layers 
for multi-shell scatter-coefficients and 6^: Nl = 370, 
layer spacing: rj = R -{- Arj^-^ with j = 1, . . . , Nl and 
Ar = R/50. Number of angles for the beam shape coeffi- 
cients: Na — 6000, mmax = 30, nmax givcu by scattnlay. 
The bottom row of the Fig. 4 shows maps of the ab- 
sorbed power Pabs- In case of a point- absorber these 
maps corresponds to the point-spread-functions (PSFs). 
These maps yield the following effective gaussian param- 
eters: cjo,^ = 281 nm, zji^d — 570 nm, Cl;o,/i = 233 nm, 
ZR,h = 468 nm. 



F. Temperatures and Cross- Sections from the 
calculations for Figure 2 

The power measured with the power-meter and con- 
verted to peak-to-peak power was Ppm = 100 /iW (see 
supplement). The power- meter aperture was the same 
as the microscope objective back-aperture, i.e. cpm = 1- 
The Gaussian GLMT calculation gives: AT§ = 96 K. 
With the expressions for c^berr 

and cg (see supplement) 
we find: Caberr = 0.384, CT,h = 0.86, cg = 1.225. The 
absorption-cross-section was found to be {zp = 0) = 
1.15 X 10^ nm^ which is close to the MiE-cross section 
{g^ = 1): crj^^e _ HQ X 10^ nm^. The exact calcu- 
lation yields AT(f = 95 K which is obtained with crjJJ^ 
and cr^g (see supplement), cr^g = 2.56 x 10~''m^, 
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1.65 X 10-^ m^. 
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1 Generalized Lorenz Mie Theory 

The Generalized Lorenz-Mie Theory (GLMT) [1,2] provides 
the total Electromagnetic fields and which solve the 
situation of the incident electric field (denoted by E^ and 
in the absence of any scatterer) upon a spherically symmet- 
ric (possibly multilayered) scattering center. The applica- 
tion of the BROMWHiCH-formalism to this situation was in- 
troduced by G. Gouesbet et al. and is here generalized for 
the purpose of finite collection- angle flux-integration. The 
field expressions, PoYNTiNG-vector expressions and far-field 
limits of the appropriate Bessel-functions involved were 
checked and adopted from [2]. For the calculation of the 
hot-particle scattering including the refractive index profile, 
the scattering coefficents Un and bn were exchanged for the 
multilayer coefficients and [8]. 

2 Integrated Fluxes 

To obtain the detected power one may calculate the in- 
tegrated PoYNTiNG-vector of the total EM-field through a 
spherical surface (see Fig. 1) at infinite distance (r oo), 
i.e. Pd=j S''dA= J S[dAvAth2{S[) = ^(^E',H';-E'^H';), 
where ^{z) denotes the real part of the complex number z. 
Also, dA = n dA with dA = r^ sin 6 dO dcj) being the surface el- 
ement in spherical coordinates and n being the surface nor- 
mal vector. The integrand may be mathematically decom- 
posed into three parts analogously to the MiE-scattering 
treatment: 

2 (Si) =di [e'^h; - E'^H';) +di[EiH; - E^Hl*) + 
^{Kh';+eih»-e;hi*-e;h';] (d 

Energy conservation dictates that -Pabs = J4n^^^'^~ ^nc + 
^sta ^ext' ^here on the right-hand side the negative sign is 
introduced to the notation to give the usual physical inter- 
pretation of the integral, i.e. Pabs being the absorbed power 
by the sphere. For this expression of energy conservation to 
hold, one always has P^'f < 0, while Pabs, P^a > 0. On the 
other hand, = 0, which means that the incident field in- 
tegral describes a flux into and out of the integration- sphere 
and a non-absorbing medium is assumed, i.e. 3 [rim) = 
where 3{) denotes taking the imaginary part. Now, let P's 
denote the integrals over the collected finite angular do- 
main, then Pd = Pine + Psca + Pext with the physically non- 
individually detectable components given by 




Figure 1: Fields and coordinates 

The interference and scatter integrals (Pext and Psca) were 
computed with Mathematica leading to the expressions 
(5) - (8). Integration over the entire polar angle of Pext.sca 
with the help of the orthogonality relations for Legendre- 
Polynomials recovers the expression for the extinction and 
scatter cross-sections ([2], eqn. (9)), i.e. Pextl6'^=7r = Pgx^ — 
-aexJo and PscAem=7i = P^d ^ <^sca/o- The intensity Iq (irra- 
diance) to use in these definitions depends on the normal- 
ization used in the description of the incident field, i.e. the 
beam shape coefficients (BSCs) used. In case of the Gaus- 
sian GLMT (section 7), Iq = 2Po/ncol refers to the focus- 
peak intensity (and the corresponding cross-sections will 
have a superscript ^, i.e. a^^^ etc.), while for the exact GLMT 
Iq refers to the on- axis intensity of the collimated gaussian 
beam illuminating the microscope objective back aperture 
(the corresponding cross-sections will have a superscript ^, 
i.e. a^. etc.). The differential cross-sections are: 
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3 EXACT BEAM SHAPE COEFFICIENTS 



O" sea (^m ) — O" sca,l (^m ) + sca,2 (^m ) (5) 

Jo n=\m=-n n'=n\^^ 

CText [dm ) = (o"ext,l (^m ) + CT ext,2 [Om )) (6) 
r dm 00 n 00 ^ 

Jo n=lm=-n n'=n\^^ ^ 

with 27r//c2 = A^ / (271^^) and ^[^| = max(l,|m|). These 
expressions have been evaluated numerically. The matrix 
elements (M/j)^ will be introduced for convenience and 
calculated for each angle Q within the sums (and prior to 
their evaluation) which helped to speed up the numerical 
calculations: 



Nn 



^§n,TE \m\ 



^§n,TM 



T\m\ 



■■(mi,i [0)y +(mi,2 [0)y 

V 2,1 Jn V 2,2 Jn 



m i-rlml ■ m M 



in(Mi,2 {e)Y+bn (mi,i {e)Y 

V 2,1 Jn V 2,2 Jn 



(7) 



(8) 



where / = and A^^ = [2n + 1] / + 1)]. The an- 
gular functions II'^' and are defined in section 4. The 
physically detectable angular power distribution (see Fig. 
2) is obtained during the numerical evaluation of the in- 
tegral in the above expressions of eqns. (5) and (6), i.e. 

Pd([6l,6l+d6l])//o=C7sca(6>+d6l) + C7ext(6>+d6l)-C7sca(6>)- 

cJext [O]- A constant incident beam offset is neglected here. 
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Figure 2: Normalized angular signal spectrum Pd[0)/ sin(^) 
for Zp = —Zr. The vertical lines correspond to the angles 
arcsin(NArf/^^) for NA^ = 0.3 and 0.75. The solid red line 
is the signal spectrum for the exact description, while the 
dashed line corresponds to an on-axis Gaussian beam cal- 
culation. 

The relation between the powers and the cross-sections 
are: P^l = a,,dn)Io and PiJ = -c7ext(7T)/o with a^A^) = 



o"sca (tt) + CTabs (tt). Explicitly, the usual integrated quantities 
are [2]: 



\n,m ^ ' 



n,m 

^abs ~ ^ext ^sca 



(9) 



3 Exact Beam Shape Coefficients 

The beam shape coefficients (BSCs) describe the incident 
fields and A similar approach as GLMT has been used 
by A. Neves et al. to model optical trapping [7]. The expan- 
sion coefficients G^^^jj.^ as given in this reference may be 
related to the BSCs {gnju/TE^ used in the framework of the 
GLMT 



m>0 
t>n,TM 

m>0 
t>n,TE 



^(^ + 1) [n-\m\)\ 
471(2^ + 1) [n + \m\)\ ^ 



(10) 



Neves' expansion coefficients describe an aberrated fo- 
cused beam. It is the resulting field when illuminating a lens 
of diameter D = Zfsmam (the microscope illumination ob- 
jective, numerical aperture NAm = ^isin(a^), effective fo- 
cal length /) by a collimated Gaussian beam of width coa . 
The over- / underfilling is then parametrized by the ratio 
D/cOa (cf. overfilling-factor Pc = {-^^ [6]). The BSCs for 
x-polarized TEM-00 focused laser beam read (for other po- 
larizations, see [7]): 



i2n-\-l)in-m) 



ikif ■n-m+l ^-imcpo 



Til 



xr 



4nn[n + 1)(^ + m)! 



i 



cos 00 


+ imS"'TM 


sin 00 




) 


sin 00 


n, 


— cos 00 





(11) 

(12) 



-- sm «2 



Im (Xa2) 



= sm «2 



P Jm [Xg^) 



where = 
i{ki cos«i 



■ +y2,0o = arctan(3//x), = -ikzzocosaz — 
- /c2Cos«2)<i, = kzposinaz = /cipsin«i 
and cos«2 = (l -(^i/^2)^sin^«iy^^. Here, Im are the 
Bessel functions of the first kind of m-th order and = 
njJ^(cos«2) and = TjJ^(cosa2) are related to the associ- 
ated Legendre polynomials (see eqn. (16)). Also the Fresnel 
transmission coefficients are introduced: 

f riz cos a2\~^ ^ friz cos«2V^ 

r^=2 1 + — -] , r^=2 — + -] (13) 

V ^1 cosai y yrii cosai J 
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4 Special Functions 



Figure 3: Integration domain of the fluxes, BSC angles 



Here, Uq, Ui and riz (see Fig. 3) are the refractive indices of 
air, immersion-oil/cover-glass and the sample medium. To 
include the intrinsic objective design aberration which can- 
cels the aberration of a glass -cover- slip under design con- 
ditions (^0 = I, Hi = Hg = ^* , dg = <i* etc.), one may in- 
troduce an additional initial aberration function in the inte- 
grand of I'^TM as a further phase i.e. an additional phase 

n, j.^ 

factor exp(-//co^*). To simplify matters, a reasonable form 
is 



^* = cos«*, ^isin«i 



(14) 



motivated by the treatment given by Woehl et al. [6] . The pa- 
rameter <i* is a fit-parameter that physically contains all oth- 
erwise non-considered aberrations that might be present 
due to non-design parameters of immersion-oil and cov- 
erglass (layer-thicknesses, refractve indices) and the addi- 
tional interface. This non- considered interface between 
coverglass and immersion- oil has low optical contrast, since 
Ug^ Hi in the experiment. This justifies those subsumma- 
tions into a single term, and eqn. (11) with I^tm including 

n, 

the factor exp(-/A:o^*) (14) yields upon conversion through 
eqn. (10) the final form of the exact incident beam expan- 
sion BSCs: 



^'TE [n + m)\ Hi 



-m ^ — im<pQ jm 



(15) 



The BSCs need to be computed numerically for positive m, 
since g'^M.ri? = (-1)"^' gTrMjE for all m. 



The angular functions appearing within the GLMT frame- 
work are: 

P^(cos0) dP^(cos^) 
n-(cos^)^ . > C(cos^)=^4- (16) 

They were derived recursively through the follow- 
ing recurrence-relations, which generalize the expressions 
given in ([8], loc. eqn. 26(a-c)) and add some trivial identities 
obtained from known identities for the associated Legendre 
polynomials : 



n« = l/|sin0| 

nx-i)"^sin"-^ 
n;;-i = cos^(2^-i)n;;:| 

2^-1 



e 



n + m — l 
cos 0U^_, U'^_, 



(17) 
(18) 
(19) 
(20) 



The fourth relation is only to be used for m < ^ — 2. For 
the derivatives of the associated Legendre polynomials, tJ^, 
one may use the following recurrence relations afterwards: 



Tl = -n{2n- l)sin 6 cos 6 n;;:| 



(21) 
(22) 
(23) 



Also the following relation for the derivatives of the 
Bessel-functions were used: 

I'm M = = [Im-l M - Im+l (X)] /2 (24) 



5 Incident and Absorbed Powers 

Per solid angle dfl = sinOdcpdO the power dP = 
I [p) pdpdcj) is contained. The Abbe sine condition relates 
the incident beam cylinder sector distance p to the an- 
gle 6 via the effective focal length / of the objective: p = 
fsinO. Therefore, dr/dO = fcos9 and the intensity is ob- 
tained from the Gaussian incident intensity profile I [p) = 
Iq exp (^-2p2/ct)^) with beam waist coa and the definition of 
an overfilling factor y = f/coa'. 

^inc= I{p)2npdpd(p 
Jo Jo 

= 271 f /2 COS ^ sin 0/o exp (-27-2 sin^^)d0 
Jo 

= Io^ [l-exp[-2fsin^em)] (25) 

here and later the subscript e denotes the reference to the 
physical back entrance aperture of the objective, and Tobj.e is 
the radius of the back aperture. The expression goes to the 
plane-wave(pw) -limit of P^^ = h^p sin^ «^ for ]^ < 1 with 
constant illumination intensity Iq filling the aperture. This 
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limit was considered in [4] . The general expression (25) may 
be used to obtain an estimate of the detected power without 
a particle, i.e. for normalization (see next section). For heat- 
ing estimates, the incident power on the particle is needed, 
and hence = arcsin(NAiii/wi) is taken, while for the de- 
tection background intensity am = Om = arcsin(NArf/^2) 
has to be used (assuming 9m < dm)- Analogously to the 
cross-sections, we may define the following two quantities: 



= arcsin(NAiii/^i) 
= arcsin(NArf/^2) 



(26) 
(27) 



obtained by considering in the exact GLMT z^-scan of crf^g 
for the aberrated case and for the un-aberrated case and 
considering the peak-ratio. Alternatively, it may be equally 
well obtained by a PSF analysis according to [6]. 

^^aberr — r , ^ — r „r.aKovr^ ^^^J 



unaberr^ 



This means, that the following absorbed powers can be 
compared and used to calculate the temperature rise in- 
duced via A To = P^\yJ (AtzkR): 



Also, cJabs = ^abs//o and thus cr^^^jG^^^ = Pabs/Cc ^hich 



may be used to calculate absorbed powers via 

T-»ill ^abs 



(28) 



In case of a power-meter detection area radius Tpm deviat- 
ing from the objective back-aperture (robj,e)> a further factor 
to be considered is: 



Gaussian: P^^^ = Ct,X, Cg Caberr 2 (7^^, (34) 



Exact: = Cr,A, CpMi^M ^ 



abs 
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(35) 



CPM--^ 



I{p)2npAp 

Jo 

l-exp(-2r^,j.,/<J 
l-exp(-2r2^/0 



I{p)2npdp 

Jo 



(29) 



Given an aperture radius Tobj.e and beam-waist C0a,e of the 
incident beam we may relate these parameters to the ef- 
fective lens parameters via Tobj = f sin 9m = /NAm/^i and 
robj,e/<Wa,e = ^obj/cOa, such that y = ^ A further factor 
of CT,h = 0.86 or CT,d = 0.83 for 532 nm (subscript h) or 635 
nm (supscript d] wavelength accounts for the transmission 
of the used objective (Olympus, UPlanSApo lOOXO). Heating 
may thus be computed via 



AtzkR' 



(30) 



To connect the exact with an approximate Gaussian cal- 
culation we note the following. The absorbed power Pabs = 
^o^abs obtained with the help of the gaussian GLMT 

absorption cross-section cr^^^. The total power i^nc^, which 
is to be used in the Gaussian approximation in order to com- 
pute the peak intensity via /q = "^PlnT^I (^<^o)' ^^ads 



-^nc^^ — (^T,Xh ^^G^^aberr^^PMPj 



,PM 
inc ' 



(31) 



wherein the factors Cg, Caberr have the following meaning: 



Cg=- 



)Obj 



TOO / /Tobj.e 

= J 7 (p) 2npdp I j 7 (p) 27rpdp (32) 



1 -exp 



^'obj,e 



1 - exp 



-2)'2NA2y 



For ^ = 1 and rii = 1.46 and NAm = 1.4 the coefficient 
takes the value Cg — 1.225. Aberrations may be included by 
introducing an intensity scaling factor. This factor may be 



A normalization is needed in order to define the relative 
photothermal signal (eqn. 40). A quantity having the units 
of an area and being proportional to the background power 
has been introduced for that purpose: crf^^. For the ex- 
act beam, af^^ from eqn. (27) and (25) should be used. 
For a Gaussian beam (on-axis) one may calculate the inte- 
grated flux of the collected beam directly, i.e. J • da with 

2S[ = [^^77^ - E'^H'*^ 12. The result may be written as 

a Cauchy-sum crf^^ = ^ ^^=1 ^fncn summands crf^^^ 
given by: 



n pB 

m=l Jo 



(36) 



[Hm^n-m+l — {—^T ^m^n-m+l] Siu^d^ 



with An='n.n-Tn and = -h . To ensure numerical 
stability for small angles, a direct recursive determination of 
An=Iln — Tn was uscd [5]: 



Ao=0, Ai = l-cos^, A2=3 + 3cos^-6cos2^ (37) 
2^-1 

=- [\ + n{n- 1) cos 61] A^_i - ; (38) 



{n-lY 



{n-lY 



The relative photothermal signal may be calculated 
quantitatively and compared to experimental values by con- 
sidering relative signals, i.e. the ratio of the time-varying 
power, which is the photothermal signal S, and the large 
constant power 7^d detected i.e. on a photodiode (voltage 
X^d): 



APpD PpD, hot ~ PpD, cold A(J ext + A(J 



at 



{O'scdi + ^^ext) lAro(P;j,Prf) (^sca + ^^ext) \l^TQ[Pyi=0,Pd) 



(39) 



(40) 
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In a modulated heating laser beam usually the effective 
heating power PpM,ef£ = Ppm/^ is detected by a power-meter 
{PpM being the peak-to-peak power used for calculations 
and static experiments), and furthermore a lock-in amplifier 
that is common to these experiments outputs a photother- 
mal (rms, root-mean-squared) voltage Igf^^ that is related to 
the peak-to-peak voltage via V^^' = V^^/ {2^/2). Therefore, 



PpM 



T/rms 



(41) 



7 Gaussian Beam Shape Coeffi- 
cients 

While the implementation of an off- axis beam according to 
[3] may be used to generate images akin to Fig. 1 of the 
article, the aberration- caused details are naturally absent 
in these calculations. Considerably easier than the off- axis 
scattering is the computation of the signal for a particle illu- 
minated on- axis by a focused Gaussian beam. In that case, 
the double and triple sums simplify and one may utilize the 
modified local approximation (MLA, [9]): 
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^abs — ^ext ^^sca- 
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